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The  oxidation  of  the  triehlorooxyphosphorus  anion  (POCI3”),  which  takes  place  in  combustion 
flames,  has  been  examined  experimentally  at  a  variety  of  temperatures  and  theoretically  via  ah  initio 
and  density  functional  methods.  The  reaction  was  examined  in  a  turbulent  ion  flow  tube  and  kinetics 
was  measured  between  300  and  626  K,  estimating  an  overall  reaction  barrier  of  1.23  keal/mol. 
Calculations  at  the  density  functional,  Moller-Plesset  second  order  perturbation,  and  coupled  cluster 
levels  of  theory  with  basis  sets  up  to  augmented  triple-f  quality  point  to  a  multistep  reaction 
mechanism  involving  an  initial  [0P(C1)3(00)J"  intermediate,  an  adduct  between  triplet  O2  with 
POCl3“,  subsequent  formation  of  a  four-membered  nonplanar  P — O — O — Cl  ring  transition  state, 
with  concomitant  breaking  of  the  P — Cl  and  O — O  bonds  to  provide  a  transient  intermediate 
[OPICOtOO*  •  -Cl]",  which,  in  turn,  converts  to  the  product  complex  (POCl2”)(CIO)  upon  formation 
of  the  Cl — O  bond  without  barrier.  The  calculated  energy  of  the  four-membered  transition  state  is 
considered  to  be  in  good  agreement  with  the  small  overall  barrier  found  by  experiment.  The  final 
step  is  responsible  for  the  large  exothermieity  of  the  reaction.  ©  2010  American  Institute  of  Physics. 
[doi:  10. 1063/1. 3299276] 


I.  INTRODUCTION 

Gcxxlings  and  his  colleagues'*^  showed  that  addition  of 
POCI3  to  flames  results  in  PO2”  and  P03“  terminal  ions.  In 
order  to  help  explain  those  results,  we  have  studied  a  variety 
of  the  ion-moleeule  steps  involved  in  the  oxidation.  Electron 
attachment  produces  POCl3“  and  POC^".^'^  The  chemistry 
of  these  ions  was  studied  with  a  number  of  simple  stable 
neutrals  and  found  to  be  unreaetive  with  species  such  as  O2, 
CO2,  and  H20.^  Therefore,  the  oxidation  from  electron  at¬ 
tachment  products  to  PO2"  and  PO3'  must  involve  minor 
species  including  unstable  neutrals.  Partial  oxidation  was 
found  in  reactions  with  ozone. ^  POCl3“,  POCU",  and  P02C1“ 
were  found  to  react  with  O3  at  greater  than  30%  of  the  col¬ 
lision  rate.  Several  of  the  products  showed  an  oxidation 
channel  and  the  reaction  of  P02Cr  with  O2  was  found  to 
produce  P03“  and  P04“.  While  the  products  in  the  latter 
reaction  arc  fully  oxidized,  the  ion  is  produced  only  in  lim¬ 
ited  abundance.  A  series  of  reactions  involving  II  atoms 
showed  that  Cl  abstraction  to  produce  HCI  was  common.^ 
That  study  also  showed  that  the  POCC  formed  in  the  POCU" 
reaction  with  H,  reacted  with  both  O2  and  O3  to  form  POj^" 
ions.  Reactions  of  O  and  N  atoms  with  PO^^CI  “  ions  also 

9  ^ 

produced  ions  with  fewer  Cl  atoms.  Recently,  a  study  in- 
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volving  PO^Cly*  with  02(0  'A^,)  showed  little  chemistry.'® 
While  these  studies  showed  pathways  to  full  oxidation,  the 
main  route  has  not  been  found. 

All  the  ion-molecule  studies  discussed  above  were  per¬ 
formed  in  a  selected  ion  flow  tube  (SIFT)."  It  is  difficult  to 
study  reactions  with  rate  constants  smaller  than 
~10“*^  cm‘^  s"'  in  a  SIFT.  Since  Ot  is  a  main  species  in 
atmospheric  plasmas,  such  as  in  the  flame  experiments,  this 
upper  limit  is  not  very  restrictive.  Slow  reactions  with  O2  can 
change  the  chemistry  due  to  its  large  abundance.  For  in¬ 
stance,  it  was  postulated  that  if  the  02*^  reaction  with  N2 

proceeded  with  a  rate  constant  of  the  order  of  10“'^  em^  s“*, 

12  13 

various  plasmas  could  be  affected.  ’  '  The  turbulent  ion  flow 
tube  (TIFT)  in  our  laboratory  is  much  more  sensitive  than 
the  SIFT  and  it  was  shown  that  the  reaction  was  much 
smaller  than  needed  to  affect  atmospheric  plasmas.'*^ 

For  the  reasons  above,  we  have  re-examined  the  reaction 
of  PC)Cl3“  with  O2  in  the  TIFT, 

P0C1,'(X-A')  +  02(X''^P 

—  POjCUlX  ‘a')  +  C10(X-II)  +  AH„,,.,.  (1) 

Reaction  (1)  is  easily  measurable  in  that  apparatus  and  fast 
enough  to  affect  atmospheric  plasmas  with  POCI3  added, 
such  as  the  experiments  from  the  Goodings  laboratory.'*^ 
Here,  we  report  kinetics  between  300  and  626  K  for  reaction 
(1),  the  upper  limit  set  by  the  stability  of  POCl3“  at  high 
temperature.  The  reaction  of  POCK"  with  O2  was  not  mea¬ 
surable.  Other  PO^Cly"  ions  are  hard  to  produce  in  the  TIFT, 
so  no  attempt  was  made  to  study  them. 
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The  detailed  study  of  the  reactivity  of  such  oxyphospho- 
rus  compounds  from  a  theoretical  standpoint  is  a  matter  that 
has  virtually  been  ignored  in  literature,  notwithstanding  the 
plethora  of  recent  experimental  data.  Early  work  of  such  type 
has  mostly  been  limited  to  the  structures  and  energies  of 
reactants  and  products^  and  other  studies  include  interaction 
energies  of  P02C12"  with  H2O  and  HCl,*"*  as  well  as  hydroly¬ 
sis  reactions  of  the  phosphate  anion.  In  our  case,  reaction 
(1)  is  unusual  in  that  the  strong  O2  bond  is  broken. 
Therefore,  in  order  to  gain  insight  into  the  reaction  mecha¬ 
nism,  detailed  theoretical  [ab  initio  and  density  functional 
theory  (DFT)]  calculations  have  also  been  performed  to 
complement  the  experimental  study. 

li.  EXPERIMENTAL  METHOD 

The  experiments  were  performed  in  a  TIFT,  which  has 
been  described  in  detail  previously.  Therefore,  only  a 
brief  description  is  given  here,  with  emphasis  on  the  experi¬ 
mental  details  relevant  to  the  present  study.  Ions  were  created 
in  a  sidcarm  with  a  corona  discharge.  The  source  region  was 
separated  from  the  flow  tube  by  a  1  mm  orifice  to  prevent 
backstreaming  into  the  source  and  unwanted  chemistry  in¬ 
volving  N  (from  the  discharge)  and  N2  from  occurring.  The 
source  gas  consisted  of  a  small  amount  of  POCI3  added  to 
Ar.  The  source  produced  mainly  a  mixture  of  POCI3"  and 
PCXri2“  by  electron  attachment.*^"^  Upon  exiting  the  ion 
source,  the  ions  were  carried  downstream  through  a  2.54  cm 
diameter  flow  tube  by  a  large  flow  of  NS  obtained  from  liq¬ 
uid  N2  boil-off.  The  ions  were  sampled  through  a  nose  cone 
and  skimmer  into  a  quadrupole  mass  spectrometer.  The  en¬ 
tire  in.strument  can  be  heated  up  to  700  K.  O2  was  added 
about  60  cm  from  the  nose  cone  and  thoroughly  mixed  in  the 
turbulent  flow  by  fan  shaped  devices  near  the  O2  inlet.  For 
these  experiments,  the  instrument  was  operated  at  100  Torr. 
This  pressure  was  about  200  times  higher  than  that  used  in 
the  SIFT,  resulting  in  longer  reaction  times  and  larger  abso¬ 
lute  number  densities  of  neutral  reagents.  This,  in  turn,  in- 
crea.scd  the  sensitivity  to  small  rate  constants.  In  these  ex¬ 
periments,  the  O2  concentration  was  about  equal  to  the  total 
buffer  pressure  in  a  typical  SIFT  experiment.  Kinetics  was 
measured  at  a  fixed  reaction  time  and  the  O2  concentration 
was  varied. 

The  temperature  of  the  flow  tube  was  changed  by  a  se¬ 
ries  of  resistance  heaters  embedded  in  insulating  blankets. 
Several  zones  were  used  for  better  temperature  uniformity. 
N2  was  preheated  as  explained  previously.  The  upper  tem¬ 
perature  limit  was  determined  by  the  fact  that  POC^"  dimin¬ 
ishes  with  temperature  and  vanished  much  above  626  K. 

Figure  1  shows  three  mass  spectra.  The  top  two  panels 
show  high  resolution  spectra  taken  at  100  with  and  with¬ 
out  O2  added.  At  this  temperature,  the  dominant  ion  in  the 
absence  of  O2  was  POCI3"  and  four  isotopomers  were 
present.  POCl2~  was  also  relatively  abundant.  These  were  the 
two  main  species  formed  by  low  energy  electron  attachment 
to  POCI3,  which  shows  that  the  electron  temperature  was 
near  thermal. In  addition,  a  small  PO2CI2”  signal  was 
present  in  the  spectrum  without  added  O2  probably  due  to  a 
5  ppm  O2  impurity  in  the  N2  buffer  or  a  small  leak.  A  few 
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FIG.  1.  High-resoluiion  ma.ss  spectra  laken  ai  100  °C  (a)  with,  (b)  wiihout 
O2  added,  and  (c)  al  350  °C  wilh  low  mass  resolution  without  O2  added. 

percent  of  Cl"  was  present.  All  other  peaks  were  much 
smaller.  The  middle  spectrum  clearly  shows  that  O2  addition 
converts  PC)Cl3"  into  PO2CI2",  leaving  the  other  peaks  essen¬ 
tially  unchanged.  The  bottom  spectrum  was  taken  at  350  ®C 
with  low  mass  resolution  without  O2  added.  At  this  tempera¬ 
ture,  the  POCl^"  signal  was  small  but  still  observable.  The 
PC)Cl3“/P02Cl2~  ratio  was  about  1,  which  shows  that  the 
rate  constant  for  reaction  (1)  increased  at  higher  tempera¬ 
tures,  More  Cr  was  also  observed. 

ill.  THEORETICAL  APPROACH 

Several  trial  geometries  of  the  (PC)Cl3,02)^  system  were 
considered  and  calculated  initially  at  the  DFT  level  of  theory, 
employing  the  hybrid  Bccke  three-parameter  Lee- Yang-Parr 
exchange  correlation  functional  (B3LYP)  and  the  6-31 
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FIG.  2.  Ploi  of  ion  intensities  vs  concentration  for  data  taken  al  300  K. 

-»-G(d)  basis  set  with  the  GAUSSIAN  03  package. Stationary 
points  found  on  the  ground  doublet  (Dq)  potential  energy 
hypersurfacc  (PES)  were  characterized  by  harmonic  fre¬ 
quency  calculations.  Connectivity  of  the  optimized  transition 
state  structures  with  the  corresponding  reactants  and  prod¬ 
ucts  was  verified  by  following  the  reaction  coordinate  back¬ 
ward  and  forward,  respectively,  using  the  intrinsic  reaction 
coordinate  methcxl.^*  Once  a  qualitative  reaction  pathway 
was  established,  we  proceeded  to  obtain  more  accurate  ge¬ 
ometries  of  all  the  structures  using  the  spin-unrestrieted  sec¬ 
ond  order  M0llcr-Plcssct  perturbation  theory  (MP2).  For  the 
onc-clcctron  ba.sis,  wc  employed  the  correlation-consistent 
basis  .sets  of  double  and  triple-^  quality,  aug-cc-pV(D  +  d)Z 
and  aug-cc-pV(T-i-d)Z.^^’^^  Diffuse  functions  (“aug")  arc 
crucial  for  correctly  describing  the  anionic  character  of  the 
system,  while  the  “-f-d"  basis  was  deemed  necessary  because 
of  the  presence  of  the  second-row  elements  P  and  Cl,  for 

which  it  has  been  previously  shown  that  the  extra  “d''  func- 

23 

tion  cures  significant  convergence  issues  as  well  as  the 
BSSE  especially  in  the  DZ  and  TZ  bases. 

To  further  increase  the  quality  of  our  findings,  single 
point  coupled  cluster  singles  and  doubles -1- a  perturbative  es¬ 
timate  of  triples  [CCSD(T)]  calculations  were  carried  out  on 
the  MP2/aug-cc-pV(T-i-d)Z  geometries  using  the  same  basis 
set.  All  correlated  calculations  keep  the  23  lowe.st-in-encrgy 
orbitals  (corresponding  to  the  inner,  non  valence  shells  of  P, 
Cl,  and  O)  as  “frozen”  (core). 

As  mentioned  already,  wc  have  used  the  GAUSSIAN  03 
software  package  along  the  theoretical  study. Some  initial 
PES  scans  were  performed  with  the  help  of  the  GRRM  eode.^^ 

IV,  RESULTS 
A.  Experimental 

Figure  2  shows  a  plot  of  ion  intensities  versus  O2  con¬ 
centration  for  the  data  taken  at  300  K.  In  this  data  set, 
PC)Cl3"  and  P(Xri2”  start  at  about  equal  intensities  and  the 
former  decreases  by  over  an  order  of  magnitude  and  the  lat¬ 
ter  is  essentially  constant.  PO2CI2"  shows  a  sharp  rise.  The 
300  K  rate  constant  is  1.25  X  10”*“^  cm^  s"*.  The  difference 


FIG.  3.  Arrhenius  plot  of  the  results  for  five  total  temperatures  over  a 
300-626  K  temperature  range. 

in  slopes  between  the  two  ions  shows  that  small  rate  con¬ 
stants  can  easily  be  measured  in  the  TIFT.  No  other  ion  in¬ 
tensity  increased,  and  therefore  the  only  product  is  PO^CU”. 
Neutral  prcxiucts  are  not  detected  and  arc  inferred  from  ther¬ 
mochemistry.  To  confirm  the  result,  we  measured  rate  con¬ 
stants  with  multiple  isotopomers  and  found  the  same  rate 
constant  within  a  few  percent.  The  error  is  estimated  to  he 
±25%. 

The  measurements  were  repeated  for  a  total  of  live  tem¬ 
peratures  over  the  temperature  range  of  300-626  K.  An 
Arrhenius  plot  of  the  results  is  shown  in  Fig.  3.  The  rate 
constant  increa.ses  slowly  with  temperature,  reaching  a  value 
of  3.73  X  10"*“*  cm^  s“‘  at  626  K.  The  estimated  relative  er¬ 
ror  is  about  15%  and  the  threefold  increase  over  the  entire 
temperature  range  is  well  outside  the  error.  An  activation 
energy  of  1.23  kcal  mol"*  is  derived  and  the  pre-exponential 
factor  is  approximately  10“*^  cm^s“*.  Reaction  (1)  is  exo¬ 
thermic  by  about  41  kcal  mol“*  (vide  infra).  Therefore,  these 
results  clearly  show  that  there  is  a  dynamical  barrier  to  the 
reaction. 


B.  Theoretical 

A  distinct  feature  of  ion-molecule  reactions  of  {A 
type  is  that  there  are  always  two  different  reactant 
(and  product)  channels,  cither  A^^“-fB  or  A-fB'*'^”;  therefore, 
the  method  chosen  should  at  least  be  able  to  predict  the  cor¬ 
rect  ordering  of  ionization  potentials  or  electron  affinities  of 
A  and  B."^  It  may  get  even  more  complicated  when  there  arc 
many  energetically  low-lying  excited  states  or  isomers  of  A 
and  B  which  could  interfere  along  the  reaction  pathway. 

A  number  of  possible  reactants  and  products  of  this  type 
have  been  considered  for  the  {POCI3-FO2— ►PO2CI2+CIO}'’ 
reaction  on  the  ground  doublet  hypersurfacc  and  the  relevant 
geometrical  structures  and  energies  at  the  B3LYP/6-31 
+  G(d)  level  arc  shown  in  Table  I  and  Fig.  4.  Detailed  dis¬ 
cussion  of  the  nature  of  each  of  these  isomers  is  beyond  the 
scope  of  this  work.  The  experimental  electron  affinities  of 
PCK2I3  and  O2  are  1.4  and  0.451  eV,  respectively,^^"*  so  the 
lowest  energy  reactants  arc  PC)Cl3"(X  "A')-f02(X 
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TABLK  I.  Total  energies  E(E^),  ZPE  (kcal/mol),  ZPE-correctcd  relative  energies  AE  (kcal/mol),  Mulliken  atomic  charges  q,  and  spin  densities  (p)  of  possible 
reactants  and  products  for  the  [POCl^+Oi]"  reaction  system  at  the  B3LYP/6-31  +G’  level  of  theory. 


E 

ZPE 

AE 

qi.(pp) 

9chi)(Po(I)) 

9o(2)(/\X2)) 

9a(i)(Pci(i)) 

9a(2)(pri(2)) 

9ci(3)(Pa(.^)) 

POCl3_So 

-1797.242  16 

5.9 

0.0 

0.36 

-0.42 

0.02 

0.02 

0.02 

POCb.To 

-1797.133  29 

4.3 

66.7 

0.58(0.48) 

-0.43(0.23) 

-0.02(0.26) 

-0.04(0.14) 

-0.09(0.88) 

POCl3‘_Do 

-1797.333  53 

4.1 

0.0 

0.46(0.30) 

-0.50(0.03) 

-0.32(0.22) 

-0.32(0.22) 

-  0.32(0.23) 

POCV_D, 

-1797.302  66 

4.0 

19.3 

0.36(0.06) 

-0.44(0.20) 

-0.29(0.04) 

-0.29(0.04) 

-0.34(0.67) 

POCI3  _D: 

-1797.301  62 

3.9 

19.8 

0.47(-0.01) 

-0.46(0.29) 

-0.33(0.03) 

-0.33(0.03) 

-0.34(0.67) 

POeV^Q, 

-1797.218  24 

3.4 

71.6 

0.30(1.20) 

-0.55(0.26) 

-0.15(0.62) 

-0.26(0.26) 

-0.34(0.67) 

PO^Clj^Do 

-1412.202  28 

6.2 

0.0 

0.78(-0.09) 

-0.38(0.52) 

-0.38(0.52) 

-0.01(0.02) 

-0.01(0.02) 

P02Cl2_D, 

-1412.171  07 

6.1 

19.5 

0.96(-0.01) 

-0.49(0.00) 

-0.49(0.00) 

0.04(1.00) 

-0.02(0.00) 

P02Cl2_D2 

-1412.144  99 

5.7 

35.5 

0.82(0.53) 

-0.37(-0.07) 

-0.51(0.23) 

0.13(0.16) 

-0.07(0.15) 

PO2CI:  D, 

-1412.108  84 

5.2 

57.6 

0.81(-0.01) 

-0.51(0.10) 

-0.40(0.01) 

0.19(0.00) 

-0.08(0.91) 

PO.Clj.Qo 

-1412.087  74 

5.0 

70.6 

0.84(0.53) 

-0.48(0.38) 

-0.31(0.83) 

0.02(1.05) 

-0.07(0.21) 

po2ar_So 

-1412.394  07 

6.8 

0.0 

0.81 

-0.62 

-0.62 

-0.28 

-0.28 

PO^Cl,"  To 

-1412.291  40 

5.3 

62.9 

0.73(0.63) 

-0.61(0.19) 

-0.46(0.39) 

-0.35(0.67) 

-0.31(0.13) 

O.-To 

-150.327  577 

2.3 

0.0 

00(1.0) 

0.0(1.0) 

02_Sn 

-150.266  052 

2.3 

38.6 

0.0 

0.0 

or_Do 

-150.349  367 

1.7 

0.0 

-0.5(0.5) 

-0.5(0.5) 

C10_Do 

-535.300  093 

1.2 

0.0 

-0.3(0.75) 

0.3(0.25) 

ClO.Qo 

-535.233  046 

0.7 

41.5 

0.1(1.96) 

-0.1(1.04) 

CIO"  So 

-535.383  379 

0.9 

0.0 

-0.8 

-0.2 

C10-_To 

-535.368  596 

0.4 

8.8 

-0.3(1.62) 

-0.7(0.38) 

Even  though  B3LYP  overestimates  the  elcetron  alfinity  of 
P(XI!l3  by  about  1  eV  (Table  I),  we  eonsider  this  acceptable 
for  our  initial  qualitative  geometry  scans,  especially  since  the 
correct  ordering  of  the  electron  affinities  is  predicted,  B3LYP 
favors  a  symmetrical  distorted-tetrahedron  POCI3"  struc¬ 


ture;  however,  at  higher  levels  of  theory,  its  symmetry  is 
reduced  to  Cj,,  but  the  effect  of  reduced  symmetry  is  not  very 
large.  The  other  structures  of  POCI3”  eorrespe^nd  to 
[POCl2*‘*Cl]"  or  [POCl  -Cl  -Cll”  complexes  and  are 
about  20  keal/mol  or  more  higher  in  energy.  The  ground  slate 


2.456 

^•^1.218 

O,T,(0.0) 

CIO  D,  (0.0) 

^1.750 

1.511^ 

po,ci;t,  (+e2.fl) 

0;  D,  (0.0) 

010*8,(0.0) 

FIG.  4.  The  optimized  structures  (bond  distances  in  angstroms  and  bond  angles  in  degrees)  and  relative  energies  (with  ZPE  correction  in  kcal/mol)  of  possible 
reactants  and  products  of  the  {P0Cl3,02}"  system  at  the  B3LYP/6-31 +G(d)  level  of  theory. 
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(•>  (b) 

FIG.  5.  Atomic  numbering  and  schematic  representation  of  the  two  possible 
approaches,  (a)  axial  and  (b)  equatorial,  of  the  On  molecule  to  the  POCI3" 
ion. 

of  neutral  POCI3  (POCl3_So)  is  characterized  by  a  polarized 
P— ►O  bond  as  Mulliken  charge  distributions  indicate 
(Tabic  I).  No  significant  O  >P  p7r-d7r  backbonding  is  ob¬ 
served.  Comparing  these  charges  with  the  ones  in  POCl^", 
we  sec  that  the  added  negative  charge  is  concentrated  on  the 
three  Cl  ligands  of  the  phosphorus  atom,  leaving  the  polarity 
of  the  P — O  bond  practically  intact,  and  especially  the  anti¬ 
bonding  nature  of  the  singly  occupied  molecular  orbital 
(SOMO) 

orbital  results  in  the  elongation  of  the  three  P — Cl  bonds  by 
about  0.3  A  in  the  anion  versus  the  neutral  (Fig.  4).  The 
corresponding  natural  bond  order  charges  show  a  similar 


situation,  but  with  even  more  polar  P — O  bonds  (qp 
=-i-l.59,  qo=-0.98  for  POCI3,  and  qp=-i-1.41,  qo=-1.04  for 
POCI3").  This  is  also  confirmed  by  atomic  spin  density  dis¬ 
tributions  (Table  I),  indicating  that  most  of  the  excess  spin  is 
located  on  the  three  Cl  atoms.  A  remaining  0.3  of  spin  re¬ 
mains  on  the  P  atom,  which  will  facilitate  the  approach  and 
bonding  of  triplet  O2  in  the  first  step  of  the  reaction  (vide 
infra).  In  both  neutral  and  anion,  the  in  situ  P  atom  finds 
itself  in  its  ground  *^5  state,  offering  its  three  lone  pair  elec¬ 
trons  to  the  three  Cl  atoms  and  donating  its  (3s-likc)  lone 
electron  pair  to  a  *D  oxygen  atom.^^  The  following  valcncc- 
bond-Lewis  scheme  attempts  to  condense  these  bonding 
characteristics  in  a  qualitative  way. 


..6 

••o 


Chemical  intuition  also  tells  that  a  doublet  POCl3~  ion 
should  react  with  triplet  (but  not  with  singlet)  Oi  by  forming 
a  pair  of  electrons  resulting  in  an  overall  doublet  state.  Per¬ 
pendicular  approach  of  the  X  oxygen  molecule  to  the 
PC)Cl3~  positively  charged  P  atom  is  expected  to  be  attractive 


TS2_«q  i26:i*  P_CompL»q  PO,Cir(8J  CIO(DJ 


FIG.  6.  Geometrical  structures  of  reactants,  intermediates,  and  products  of  the  POCIT+Oi  reaction  in  the  equatorial  reaction  pathway  at  the  B3LYP/6-31 
+G(d)  (upper  numbers)  and  MP2/aug-cc-pV(T+d)Z  (lower  numbers)  levels  of  theory. 
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pocC(Dj  +  o/rj 


pooa=ioo.r 

105,3* 


LM1  ui 


I 

% 

R_Compl_«x 


LMO  IX 


FIG.  7.  Getimetrical  structures  of  reactants,  intermediates,  and  products  of  the  POCl^  +0^  reaction  in  the  axial  reaction  pathway  al  the  B3I.YP/6-3I 
+G(d)  (upper  numbers)  and  MP2/aug-cc-pV(T+d)Z  (low-er  numbers)  levels  of  Iheory, 


due  an  ion-dipole-indueed  Coulomb  attraction,  gradually 
leading  to  a  o*  electron  coupling  between  a  7r*0  electron  and 
the  lone  electron  of  POCI3”  according  to  this  scheme. 


O" 

1  Cl 


[o— P(a)^o-6] 


As  shown  in  Fig.  5,  two  different  orientations  were  con¬ 
sidered  for  the  attack  of  the  triplet  oxygen  molecule  to  the 
doublet  POCI3”  ion  in  the  overall  doublet  spin  state,  resulting 
in  axial  (ax)  and  equatorial  (eq)  and  geometry  configura¬ 
tions.  In  the  former,  O2  approaches  parallel  to  the  P — O  axis 
of  POCI3’,  while  in  the  latter,  it  does  so  in  a  parallel  fashion 
to  one  of  the  P— Cl  axes.  For  both  orientations,  similar 
pathways  were  determined  (Figs.  6-8)  involving  the  forma¬ 
tion  of  at  least  two  transition  state  complexes  and  two  inter¬ 
mediates  in  addition  to  van  der  Waals  reactant  and  product 


complexes.  The  results  for  the  equatorial  and  axial  orienta¬ 
tions  at  different  levels  of  theory  are  tabulated  in  Tables  II 
and  III  and  are  depicted  in  Figs.  6  and  7,  respectively.  Tables 
IV  and  V  present  Mulliken  atomic  charge  and  spin  density 
distributions  at  the  B3LYP/6-31  -fG*  level  along  both  reae- 

I#  #4- 

■<j^ 

LMO_«q_HOMO  T81_«r.HOIIO  LM1_«q.HOMO 

4^  A#® 

Tt2.«R.HOMO  P,Contpl_«q_H(HiO 

FIG.  8.  Contour  plots  (thre.shold=().05)  of  HOMO  for  selected  chemical 
species  of  interest. 
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TABLE  II.  Stationary  point  relative  energies  (kcal/mol)  at  ditTerent  levels  of  theory  for  the  equatorial  reaction 
pathway  without  and  with  zero- point  energy  differences  (AZPE). 


B3LYP/6-31 +G*  UMP2/aug-cc-pV(T+d)Z  UCCSD(T)/aug-cc-pV(T+d)Z“ 


AE 

AE+AZPE 

AE 

iE+AZEE" 

AE 

ae+azpe' 

Reactants 

0.0" 

0.0 

O.O'^ 

0.0 

0.0" 

0.0 

R  Compl  _eq 

-0.5 

-0.4 

-2.2 

-1.2 

-2.0 

-1.0 

TS0_eq 

7. 1 

8.2 

LM0_eq 

3.3 

5.5 

-1 2.2 

-9.2 

-14.3 

-11.2 

TSl  eq 

1 2.2 

1 3.8 

1.6 

4.1 

-3.5 

-0.9 

LM I  _eq 

0.5 

2.7 

-14.6 

-6.9 

-14.0 

-6.3 

TS2_eq 

2.6 

4. 1 

-5.9 

-3.3 

-11.6 

-9.0 

P_Compl  eq 

-30.4 

-28.3 

-50.6 

-48.1 

-53.7 

-51.2 

Products 

-20.7 

-I9.2 

-39.4 

-37.6 

-43.1 

-41.3 

"At  UMP2/aug-cc-pV(T+d)Z  geometry. 
*^ith  UMP2/aug-cc-pVfD+d)Z  frequencies. 
"H=- 1947.661  105  E^. 

‘^£=-1945.457  112  E^. 

’^£=-1945.575  327  E^ 


tion  pathways.  Depictions  of  the  highest  (singly)  occupied 
molecular  orbitals  (HOMOs)  for  certain  species  of  interest 
arc  presented  in  Fig.  8.  Finally,  Fig.  9  shows  the  schematic 
energetics  for  both  reaction  pathways. 

V.  DISCUSSION 

As  expected  and  shown  in  scheme  2,  a  weak  ion- 
moleculc  van  dcr  Waals  reactant  complex  (R_Compl)  is  con¬ 
verted  to  an  side-on  [0P(C1)3(00)]”  intermediate  (LMO)  for 
both  equatorial  and  axial  orientations.  This  step  bears  many 
similarities  to  the  known  reaction  mechanism  between  alkyl 
cations  and  As  mentioned  already,  LMO  is  formed 

when  02(X  points  one  of  its  two  tt*  unpaired  electrons 

toward  the  positively  charged  P  atom  and  eventually  forms  a 
single  if  bond  as  scheme  2  shows.  The  formation  is  barrier¬ 
less  at  high  levels  of  theory,  although  a  transition  state  (TSO) 
was  found  at  the  B3LYP/6-31 +G(d)  level  which  could  not 
be  found  at  the  MP2  level.  According  to  Mulliken  charge  and 
atomic  spin  density  distributions  (Tables  IV  and  V),  an  elec¬ 
tron  flux  from  the  chlorine  atoms  toward  dioxygen  via  phos¬ 


phorus  takes  place.  Indeed,  the  three  chlorine  atoms  lose 
about  ^0.35e”  upon  formation  of  the  P — O — O  bond.  The 
P  atom  also  loses  0.1  Se*",  enhancing  the  pK^arity  of  the  pre¬ 
existing  P — O  bond.  The  reactants  and  the  reactant  complex 
(R_Compl)  have  a  total  of  three  unpaired  electrons  coupled 
into  a  doublet  {aap)^  which  upon  formation  of  LMO  are 
recoupled  into  a  single-electron  doublet,  with  the  spin  lo¬ 
cated  on  the  O — O  moiety,  essentially  on  the  other  tt*  or¬ 
bital  (Tables  IV  and  V)  and  formally  vertical  to  the 
P — O — O  plane. 

When  zero-point  energy  differences  (AZPES)  are  taken 
into  account,  the  formation  of  this  P — O  bond  is  exothermic 
by  11.2  kcal/mol  for  LM0_eq,  but  only  0.9  keal/mol  for 
LM0_ax  at  the  CCSD(T)/aug-ce-pV(T+d)Z  level  of  theory. 
There  is  also  an  *^0.01  A  difference  in  the  length  of  the 
formed  P — O  bond  between  these  two  isomers.  Thus,  an 
axial  effect  seems  to  be  significant.  Mulliken  charges  reveal 
that  the  O — O  moiety  attracts  less  electron  charge  in  the 
more  weakly  bound  LM0_ax  (0.26e“  versus  0.4e'”).  The  ef¬ 
fect  of  the  axial  ligand  is  demonstrated  by  its  Mulliken 


TABLE  111.  Stationary  point  relative  energies  (kcal/mol)  at  different  levels  of  theory  for  the  axial  reaction 
pathway  without  and  with  zero-point  energy  differences  (AZPE). 


B3LYP/6-31 +G*  UMP2/aug-cc-pV(T+d)Z  UCCSD(T)/aug-cc-pV(T+d)Z" 


AE 

AE+AZPE 

AE 

AE+AZPE^ 

AE 

ae+azpe’ 

Reactants 

0.0" 

0.0 

0.0*^ 

0.0 

0.0" 

0.0 

R_Compl_ax 

-0.3 

1 

p 

-2.2 

-1.2 

-1.9 

-0.9 

TSO  a.\ 

9.4 

10.9 

LM0_ax 

9.4 

11.2 

-2.3 

0.8 

-4.0 

-0.9 

TS 1  _ax 

12.0 

13.7 

2.1 

4.5 

-3.4 

-1.0 

LMI  ax 

1.6 

3.8 

-12.6 

-6.2 

-12.3 

-5.9 

TS2_ax 

3.1 

4.5 

-4.7 

-2.2 

-10.5 

-8.0 

P_CompLax 

-30.4 

-28.3 

-50.6 

-48.1 

-53.7 

-51.3 

Products 

-20.7 

-19.2 

-39.4 

-37.6 

-43.1 

-41.3 

"At  UMP2/aug-cc-pV(T+d)Z  geometry. 
*^Vith  UMP2/aug-cc-pV(D+d)Z  frequencies. 
"£=-1947.661  105  E^. 

‘‘E=-I945.4S7  112  E,. 

"E=-I94.‘>.575  727  E*. 
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TABLE  IV.  Mulliken  alomic  charge  and  spin  densiiy  disiribulions  (in  parentheses)  al  ihe  B3LYP/6-31 +G* 
level  of  theory  {equatorial  reaction  pathway). 


Atom 

R_compl_eq 

LM0_eq 

TSI.eq 

LMLeq 

TS2_eq 

P_compl_eq 

I 

P 

0.42(-0.29) 

0.57(-0.01) 

0.62(0.00) 

0.79(0.00) 

0.80(0.00) 

0.78(0.00) 

2 

0 

-0.48(-0.03) 

-0.55(0.00) 

-0.50(0.00) 

-0.58(0.00) 

-0.59(0.00) 

-0.62(0.00) 

3 

Cl 

-0.32(-0.22) 

0 

0 

0 

1 

-0.27(0.00) 

-0.16(0.00) 

-0.18(0.00) 

-0.25(0.00) 

4 

Cl 

-0.31  (-0.22) 

-0.16(-0.01) 

-0.51(0.18) 

-0.52(0.51) 

-0.43(0.43) 

0.34(0.28) 

5 

Cl 

-0.32(-0.23) 

-0.29(0.00) 

-0.05(0.01) 

-0.07(0.00) 

-0.11(0.00) 

-0.24(0.00) 

6 

0 

0.03(1.03) 

-0.15(0.35) 

-0.22(0.13) 

-0.22(0.04) 

-0.30(0.07) 

-0.62(0.00) 

7 

0 

-0.02(0.97) 

-0.25(0.67) 

-0.07(0.67) 

-0.24(0.45) 

-0.18(0.49) 

-0.39(0.71) 

charge;  in  both  axial  and  equatorial  cases,  the  axial  ligand 
(Cl  or  O)  becomes  more  negatively  eharged  due  to  the  attack 
of  Ot  to  P. 

The  reaetion  from  LMO  proeeeds  similarly  for  both  axial 
and  equatorial  orientations.  A  four-membered  nonplanar 
P*“0 — 0---C1  ring  transition  state  (TSl)  is  formed,  essen¬ 
tially  upon  the  transfer  of  one  Cl  atom  from  P  to  O,  and  the 
subsequent  breaking  of  the  P — Cl  bond,  to  result  in  the  for¬ 
mation  of  the  .second  intermediate  complex 
{  OP(Cl)200*  • -Cl]  "  (LMl)  with  a  partially  formed 
— 0“-C?  bond  (Figs.  6  and  7).  At  the  B3LYP/6-31 +G(d) 
level,  TSl  energetically  lies  significantly  above  the  reactants 
(~14  keal/mol)  for  both  pathways.  At  the 
MP2/aug-cc-pV{D+d)Z  level,  both  TSl_eq  and  TSl_ax 
have  one  imaginary  frequeney  eaeh  (114/  em"‘  and 
174/  cm~\  respectively)  with  the  much  reduced  energy  of 
~4-5  kcal/mol  above  the  reactants.  The  energies  of  both 
TSl_eq  and  TSl_ax  beeome  slightly  negative,  from  —1.0  to 
—0.9  kcal/mol  relative  to  the  reactants  at  our  best 
CCSD(T)/aug-cc-pV(T-i-d)Z//MP2  level.  These  are  the 
highest  points  on  the  doublet  PES  and  eorrespond  to  the 
overall  barriers  of  the  reaction.  There  seems  to  be  no  prefer¬ 
ence  between  the  equatorial  and  the  axial  pathways.  This 
barrier  ean  be  eompared  with  the  experimentally  observed 
barrier  for  the  reaction  of  1.23  kcal/mol.  Although  the  sign  is 
wrong,  the  very  small  calculated  barrier  is  taken  to  be  in 
satisfactory  agreement  with  the  experiment.  The  error  of 
~2.2  keal/mol  between  the  ealculated  and  the  experimental 
barriers  may  be  attributed  to  the  use  of  MP2  geometry,  MP2 
frequencies,  and  the  basis  set  ineompleteness  in  the 
CCSD(T)  calculations. 

The  final  step  of  the  reaction  takes  place  from  LMl  via 
TS2,  where  the  Cl  atom  with  a  long  O*  •  -  Cl  bond  is  captured 
by  the  terminal  oxygen  atom  and  the  elongated  and  weak¬ 


ened  O — O  bond  breaks  forming  an  ion-molecule  complex 
P_Compl  between  the  products:  C10(X  ^11)  and  P02CL~. 
Mulliken  atomic  charges  and  spin  densities  for  both  equato¬ 
rial  and  axial  orientations  confirm  these  observations  (Tables 
IV  and  V).  MP2/aug-ee-pV(D+d)Z  gives  a  large  imaginary 
frequency  of  932/  em"‘  for  TS2_eq  and  898/  cm"'  for 
TS2_ax,  indicating  that  TS2  is  a  real  transition  state  at  these 
levels,  with  the  energy  requirement  of  3-4  keal/mol  from 
LMl  and  is  energetically  below  the  reactants.  However,  at 
the  CCSD(T)/aug-ee-pV(T+d)Z  MP2evel,  when  AZPE  is 
taken  into  account,  TS2  becomes  lower  in  energy  than  LMl 
and  the  transition  state  vanishes.  The  change  in  energy  from 
LMl  to  the  product  complex  P_Compl  is  very  negative, 
dominating  the  large  overall  exothermicity  of  the  reaction. 
The  product  complex  P_Compl  dissociates  into  the  isolated 
products  POCL^+CIO  with  a  energy  requirement  of 
--10  keal/mol  and  an  overall  exothermicity  of  41.3  kcal/ 
mol  at  the  CCSD(T)  level  (Tables  III  and  IV). 

VI.  SUMMARY 

The  POCI^'+Ot  reaction  has  been  studied  from  both 
experimental  and  theoretical  viewpoints.  The  main  products 
of  the  reaction  are  P02Cl2"(Y  'A')  +  C10(X  “II).  The  reaction 
was  examined  in  a  TIFT,  and  kinetics  was  measured  between 
300  and  626  K.  These  measurements  gave  an  overall  reaction 
barrier  of  1.23  kcal/mol. 

Examination  of  the  reaction  pathways  at  different  levels 
of  theory  points  out  to  a  multistep  reaction  mechanism,  in¬ 
volving  an  initial  [0P(C1)3(00)]“  intermediate,  an  adduet 
between  O2  with  POCl^”,  subsequent  formation  of  a  four- 
membered  nonplanar  P — O — O — Cl  ring  transition  state, 
and  with  concomitant  breaking  of  the  P — Cl  and 
O — O  bonds  to  provide  a  transient  intermediate 


TABLE  V.  Mulliken  alomic  charge  and  spin  density  disiribulions  (in  parenthe.ses)  at  the  B3LYP/6-31 +G*  level 
of  thei^ry  {axial  read  ion  pathway). 


Alom 

R  com  pi  ax 

LMO  ax 

TSl  ax 

LMl  ax 

TS2_ax 

P  compl  ax 

1 

P 

0.44(-0.3l) 

0.68(-0.06) 

0.62(0.00) 

0.79(0.00) 

0.79(0.00) 

0.78(0.(X)) 

2 

0 

-0.51(-0.03) 

-0.68(0.t)0) 

-0.49(0.01) 

-0.54(0.00) 

-0.56(0.00) 

-0.62(0.00) 

3 

Cl 

-0.32(-0.22) 

-0.25(-0.04) 

-0.28(0.00) 

-0.16(0.00) 

-0.18(0.01) 

-0.25(0.00) 

4 

Cl 

-0.31  (-0.22) 

-0.24(-0.03) 

-0.50(0.19) 

-0.51(0.53) 

-0.43(0.45) 

0.34(0.28) 

5 

Cl 

-0.31  (-0.22) 

-0.24(-0.03) 

-0.05(0.00) 

-0.13(0.00) 

-0.14(0.00) 

-0.24(0.00) 

6 

0 

0.04(1.03) 

-0.06(0.49) 

-0.17(0.14) 

-0.21(0.03) 

-0.29(0.07) 

-0.63(0.00) 

7 

0 

-0.03(0.97) 

-0.20(0.68) 

-0.13(0.66) 

-0.24(0.44) 

-0.19(0.47) 

-0.39(0.71) 
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FIG.  9.  Comparison  of  ener|5etics  (with  ZPE  correction)  for  the  equatorial 
(.straight  lines)  and  axial  (wiggly  lines)  pathways  at  different 
levels  of  theory  (red; B3LYP/6-31G*,  green :MP2/aug-cc-pV(T+d)Z, 
hlack:CCSD(T)/aug-cC“pV(T+d)Z  at  MP2  geometry).  The  da.shed  line 
shows  the  experimental  value  of  the  energy  barrier  (1.23  kcal/mol). 

[0P(CI)200*  •  Cl]"*,  which  converts  to  the  product  complex 
(POCl2“)(C10)  upon  formation  of  the  Cl — O  bond  without 
barrier.  The  calculated  energy  of  the  four-membered  transi¬ 
tion  state,  the  highest  point  on  the  potential  profile  at  the 
highest  level,  is  ~-I.O  kcal/mol,  which  is  considered  to  be 
in  good  agreement  with  the  small  overall  barrier  by  experi¬ 
ment.  The  final  step  is  responsible  for  the  large  cxothermicity 
of  the  reaction. 
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